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ABSTRACT 

Recent theories suggest planetesimal formation via streaming and/or gravitational in- 
stabilities may be triggered by localized enhancements in the dust-to-gas ratio, and 
one hypothesis is that sufficient enhancements may be produced in the pile-up of small 
solid particles inspiralling under aerodynamic drag from the large mass reservoir in 
the outer disc. Studies of particle pile-up in static gas discs have provided partial 
support for this hypothesis. Here, we study the radial and temporal evolution of the 
dust-to-gas ratio in turbulent discs, that evolve under the action of viscosity and pho- 
toevaporation. We find that particle pile-ups do not generically occur within evolving 
discs, particularly if the introduction of large grains is restricted to the inner, dense 
regions of a disc. Instead, radial drift results in depletion of solids from the outer disc, 
while the inner disc maintains a dust-to-gas ratio that is within a factor of ^ 2 of the 
initial value. We attribute this result to the short time-scales for turbulent diffusion 
and radial advection (with the mean gas flow) in the inner disc. We show that the 
qualitative evolution of the dust-to-gas ratio depends only weakly upon the parame- 
ters of the disc model (the disc mass, size, viscosity, and value of the Schmidt number), 
and discuss the implications for planetesimal formation via collective instabilities. Our 
results suggest that in discs where there is a significant level of midplane turbulence 
and accretion, planetesimal formation would need to be possible in the absence of 
large-scale enhancements. Instead, trapping and concentration of particles within lo- 
cal turbulent structures may be required as a first stage of planetesimal formation. 

Key words: accretion, accretion discs — planets and satellites: formation — proto- 
planetary discs — stars: prc-main-sequence — stars: variables: T Tauri, Herbig Ae/Be 



1 INTRODUCTION 



Uncovering the pathway by which sub-micron dust grains 
grow into gravitationally bound bodies suitable for building 
planets is a long-standing problem. Pairwise collisions be- 
tween dust particles in protoplanet ary discs are thought to 
drive rapid growth up to ~mm sizes (|Dullemond fc DominikI 
[2OO5I). Beyond this size scale, continued growth up to 
the km-scale of planetesimals is potentially frustrated by 
two distinct physical barriers. First, the measured ma- 
terial properties of solid aggregates, which favor sticking 
upon collision for small sizes, switch to favour bouncing 
l|Giittler et al.l I2OI0I : IZsom et al]|2010t ). or coUisional frag- 
mentation due to high, turbulence-driven velo city disper- 
sions (e.g. iBirnstiel. DuUemond fc Braueiil2010l ) for Sd < m. 
Even if these barriers could be surmounted, however, di- 



rect coUisi onal growth faces a s econd obstacle from head- 
wind drag (|WeidenschillineJI 19771 ) . which becomes extremely 
strong for particles in the cm to m-size range. The loss of 
substantial quantities of solids into the star is not in itself 
deleterious - indeed the relativ ely small solid mass in the 
Minimum Mass Solar Nebula (|WeidenschiUing| [l977bl ) in- 
dicates loss probably occurred - but the residence time of 
critically coupled particles is so short that it is hard to see 
how catastrophic depletion could be averted if bodies grow 
steadily through m-scales. 

These considerations have led to renewed interest in 
collective mechanisms for planetesimal formation, which ap- 
peal (typic ally) to gravita t ional instability in a den se par- 
ticle layer (|Safronovl [iQT^ : iGoldreich fc Wardlll97J ) to ef- 
fect a dynamical time scale jump from mm or cm-sizes to 
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km-scales (for a review see IChiang fc Youdinll2010l ). Sev- 
eral precursor mechanisms have been proposed that might 
act to create conditions locally favourable fo r gravitational 



insta bility. They include vertical settling (lYoudin fc Shu 



2002 ) , turbulent concentr ation (jCuzzi. Hogan fc Sharif] 
200S: Ijohansen et alj 120071 ) , and streaming instabilities 



llJohansen fc YoudinI 2007 : Ijohansen. Youdin fc Mac Lowl 
l2009bl : iBai fc Stonj bOlOal ). There is some evidence that 
these mechanisms could form surprisingly large first- 
generation planetesimals ( jj 100 km), of a size consistent 
with that derived — by combining the measured size distri- 
bution with the results of coUisio nal-grinding simulatio ns — 
for the primordial asteroid belt (JMorbidelli et al.ll2009l ). 

Which, if any, of these mechanisms is responsible for 
planetesimal formation is hard to say, since their feasibil- 
ity depends in part on the uncerta in properties of the in- 
trinsic turbulence within the disc (|Armitagd l201ll ). Some 
theoretical constraints are nonetheless possible, since most 
of the formation-via-coUapse models require local enhance- 
ment of the height-averaged dust-to-gas ratio (the local 
disc metallicity) above solar values in order to exceed the 
threshold for gravitational instability q The magnitude 
of the required enhancement varies both with the model, 
and with the properties of the gas disc, but typical values 
quoted vary between roughly 2 and 10 times standard metal- 
licity jYoudin fc Shul |2002(: Ijohansen. Youdin fc Mac Lowl 



l2009bl : iBai fc Stond 1201081 ). Eithe r a local build-up of 
solids or a local depletion of gas (|Throop fc Ballvl l2005l : 
lAlexander fc Armitagg|2007l ) would be equally effective. 

Further clues come from observations. To match Solar 
System constraints, conditions for planetesimal formation 
must extend over a significant radial range (because plan- 
etesimals formed with substantially distinct compositions) 
and persist for a large fraction of the disc lifetime. The evi- 
dence for an extended era of planetesimal formation comes 
from age dating of primitive meteorites, and in particular 
from the finding that the chondrules wi thin those meteorite s 
differ in age by as much as 2.5 Myr JAmelin et al.ll2002h . 
Furthermore, recent age-dating of iron-core meteorites sug- 
gests the era of planetesimal formation extended nearly as 
far back as the beginning of Solar-system history, as de- 
fined by the formation of calcium-aluminium-r ich inclusions 
(CAIs) JKleine et al.ll2005l : lBottke et al.ll2006l ). 

The replenishment of cold dust in debris discs (|Wvattl 

[2003) implies that planetesimal populations around other 
stars are not solely confined to the inner disc, but there 
is no observational requirement that planetesimal forma- 
tion spans the entire range of radii where gas was once 
present. Given this, one obvious way to create an enhance- 
ment of the dust-to-gas ratio at small to intermediate disc 
radii is to appeal to the large reservoir of mass typically 
present further out. If small particles spiralling inward pile- 
up as they migrate radially, then the "problem" of headwind 



^ The ICuzzi. Hogan fc Sharifil 1I2OO8I ') model of planetesimal for- 
mation may be an exception as it considers only height-local 
dust/gas enhancement, possibly achievable by vertical settling 
alone. 



drag becomes instead a necessary ingredient for planetesi- 
mal formation. In a static (inviscid) disc model, headwind 
drag acting on an initially uniform dust-to-gas ratio does 
in fact result in just such a pile-up, with peak values of 
the enhancement exceeding those needed to produce gravi- 
tatio n al collapse via vertical s edimentation (|Youdin fc Shul 
l2002l : I Youdin fc Chiand I2OO J ). Although not the focus of 
the original work, radial pile-ups would also result in highly 
favorable conditions for triggering streaming instabilities. 

Inviscid-disc models, which assume that the gas has nei- 
ther a net radial fiow nor local turbulent motions, include 
only the minimal set of potential aerodynamic effects that 
can modify the radial particle distribution. Such models also 
provide a poor representation of the g lobal disc environment , 
which is expected (and observed, e.g. lHartmann et al.lll998l ) 
to change substantially over the time period of interest for 
planetesimal formation. Our goal in this paper is to quantify 
the variation in the dust-to-gas ratio if the gas disc, instead 
of being static, evolves with time under the action of local 
turbulent transport of angular momentum and photoevap- 
oration. Headwind drag will still occur in an evolving disc, 
but the propensity of particles to pile-up will additionally 
be affected by several new processes, 

• The global evolution of the gas disc, which thins with 
time and whose outer edge must expand radially to conserve 
angular momentum as gas is accreted on to the star. 

• Radial advection with the flow of disc gas, which may 
counter the slow-down of drift that leads to pileup in the 
inner disc. 

• Turbulent diffusion. 

• Disc-clearing photoevaporation. 

The first three of these p rocesses (also explored in, e.g., 
IStepinski fc Valagea s 1996) are ultimately byproducts of the 
same turbulence within the disc, so the inclusion of evolution 
introduces fewer independent parameters into a calculation 
of particle transport than one might at first fear. Indeed, the 
primary conclusion of our work is that a range of disc models 
all produce roughly the same result: as the disc evolves, the 
outer disc becomes depleted of solids, but the lost material 
does not pile-up significantly in the inner disc. Instead, the 
inner disc maintains roughly its primordial dust-to-gas ratio 
as the disc drains away. As time progresses, the thinning gas 
supports successively smaller particle sizes across the range 
of Solar System radii. 

The plan of this paper is as follows. In !j2]we outline the 
one dimensional disc evolution and particle transport mod- 
els that we employ. These models immediately specify the 
radial dependence of several critical time scales relevant to 
particle transport, and in |31we discuss the elementary con- 
sequences of those time scales for dust evolution. In !j4l we 
show the dust distributions resulting from our simulations, 
which we compare in iJS] to the dust-to-gas enhancements 
thought necessary to lead to precipitation (form ation via 
vertical sedimentation as in lYoudin fc Shul (|2002l )) of large 
bodies, and, briefly in ^^ to observations of dust in discs 
around other stars. We discuss the implications of our re- 
sults for planetesimal formation models in ijTJ 
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2 METHODS 

2.1 Gas disc evolution 

Our simulations assume a ID, viscously evolving disc model, 
whose f = surface density profile is Eg ex R~^ at small 
radii, with a characteristic exponential fall-off from R = Rd, 




where Eg is the surface density of the disc gas, Mo,o is 
the initial disc mass, R is radial distance from the central 
star, and Rin = 0.1 AU is the inner boundary of the disc- 
model grid, fn the fiducial model, Mo,o ~ O.OSMq and Rd = 
20 AU, but simulations are also run in discs with Md,o ~ 
O.OIM0 and O.O9M0 01 Rd = 5 AU and 40 AU. 

Viscous evolu tion of the dis c follows the thin-disc model 
of disc evolution (|Pringlelll981^ on a ID grid of 600 points 
spa ced logarithmically in R fro m 0.1 to 15,000 AU. We use 
the IShakura fc SunvaevI (| 19731 ) alpha-parameterization for 
the disc viscosity. 



v (R) = acsHg = aQ,KHg 



(2) 



where v is the viscosity, Ca is the local gas sound speed, 
f^K is the local Keplerian angular velocity, and Hg = 
fl^^ \/ ksT / fj,mH , is the vertically- isothermal disc-gas scale- 
height, fcs is Boltzman's constant, and hitih is the average 
mass of a gas particle, set equal to 2.34 hydrogen atoms. In 
the fiducial model, a = 10~ , but simulations are also run 
in discs using a = 5 x 10~^ and 2 x 10~^. The disc-evolution 
model also includes photoevaporation at the end of the disc 
lifetime due to lO'*^ photo ns s~^ of EUV radiation, as in 
iHughes fc Armitag3 (|2010h . Equation (1). 

We assume that the vi scosity in Equation ^ der ives 
from a local turbulent stress (JBalbus fc Papaloizoull 19991 ) . In 
this case, both the evolution of the gas disc profile, and the 
magnitude of the particle diffusivity, Dp, share a common 
scaling with the midplane disc temperature. Dp oc v cc T. 
Due to gas infall and accretional heating, disc tempera- 
tures at early time s are much hotter than those of later 
discs (jCassenluOOlf ). To evaluate the evolution of the mid- 
plane temperature, we assume that local heating by viscous 
dissipation and stellar irradiation is balanced by radiative 
cooling. We compute these terms using standard methods 
(detailed in Appendix |X]| , assuming a constant dust-to-gas 
ratio (i.e., ignoring coupling to the evolving dust distribu- 
tion). For computational economy, when running our dust- 
transport simulations, we model the gas-disc evolution using 
a fit to the computed energy-balanced temperature evolu- 
tion. At a given time. 



T {R) = max 



Tclo 



Tau (t) 



R 

lAU 



9T(t) 



(3) 



where Tdoud = 10 K is the background temperature of the 
star-forming cloud, and Tau and qx follow an exponential 
decay in time from an initially hotter, steeper distribution 



(Tau ~ 500 K, qr ~ —0.6) to a cooler, shallower distribution 
appropriate for a passively-heated, flared disc (Tau ~ 280 K, 
qr ~ —0.5). To compute the stellar irradiation, which dom- 
inates at late times and large radii, we assume a stellar lu- 
minosity of Li, — 5I/soiar and a disc-starlight intersection 
angle set to cfi — 0.05 everywhere. We have included the 
time-evolving nature of the disc-temperature distribution 
for completeness in considering the fundamental differences 
between an evolving versus a static disc-mass distribution. 
However, at the time-scales / masses and accretion rates 
considered in our models, disc surface-density evolution is 
nearly identical using either a static or evolving tempera- 
ture distribution. We have also verifled that using a fit to 
the evolving temperature distribution, rather than calculat- 
ing it in each dust-transport simulation at every {R, t), does 
not affect our dust-distribution results. 

The disc-gas radial velocity, Vr,g , is calculated using the 
ID accretional velocity derived due to the viscous evolution 
of the disc surface-density profile. 



Egi?i/2 dR 



{u-LgR^'^ 



(4) 



This accretion-fiow velocity is predominantly inward 
but includes outward fiow at the outer edge of the expanding 
disc. The inclusion of the accretion velocity of the gas marks 
one of the primary differences bet ween our evo lving-disc 
model and the inviscid-disc model of lYoudin fc Sh u (200^. 
It is important because, for a typical gas surface density pro- 
file (E oc R~^), \vr,g\ does not decrease toward smaller disc 
radii, unlike the headwind-drag velocity. We note, however, 
that it is possible for discs to evolve without there being a 
significant accretion fiow at the d isc midplane. This situa- 
tion occurs in dead-zone models ()Gammidll996l : lArmitagd 
l2011f ) , in which vertical variation in the efficiency of angular 
momentum transport results in gas accretion that occurs 
primarily along the disc surface. Qualitatively, we expect 
that particle transport within dead-zone disc models would 
differ substantially depending on whether the scale height of 
the particle disc was large enough for the particles to expe- 
rience the (rapid) surface layer fiow. However, more detailed 
study of these models is left to others. 

The gas azimuthal velocity, w^.g, effecting headwind- 
drag on the simulation grains, is calculated at the disc 
midplane assuming pre ssure balance within the vertically- 
isothermal disc model (JTakeuchi fc Linll2002l ). For the nu- 
merical calculation, we use 



nt = 



"0,3 






Hi 

i?2 



(7E + -qr 



(5) 



where Qg is the gas angular velocity, qr is the power-law 
index of the temperature distribution used in Equation (|3]), 
and qs is the locally-calculated power-law index of the disc 
surface-density profile. 
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2.2 Particle Transport Setup 

2.2.1 Particle Motions 

We use the pa r ticle-t ransport code outlined in 
iHughes fc Armitage! (|201(J ) to evolve the dust distri- 
bution, which is modeled as an ensemble of particles. 
Each particle is subject to aerodynamic advection and 
random-walk turbulent diffusion within the disc. Particle 
trajectories are integrated following 



D^=Da 



rd {t + At) = Td (t) + {Vsrd ± Uturb) At , 



(6) 



where Vd is a particle's radial position, At the time-step, v^rd 
is a particle's radial velocity due to gas-drag and advection 
within the mean gas flow, and Uturb is the mean velocity for 
effecting random- walk turbulent diffusion of the particle en- 
semble. The random walk is asymmetric, with the weighting 
between the inward and outward steps calculated from the 
local di sc conditions so as to recov er the evolution of a fluid 
model (|Hughes fc Armitage|[2O10l 'l. 

Gas drag is calculated assuming that particles experi- 
ence Epstein drag appropriate to their size. The Epstein 
drag regime is appropriate for the vast majority of the 
disc's radial and temporal extent, though the largest grains 
considered (cm-sized) technically fall within the Stokes- 
drag regime in the innermost disc ( ^ 1 AU) at very 
early times. The radial velocity of a grain is calculated at 
each radial grid-point by iteratively solving for the steady 
velocity of particles that have attained termin al velocity 
iTakeuchi fc Linll2002l : iHughes fc Armitagellioioh . 



'^srd 



f^srd 



+ 



vl) 



CRRpgVtho 



--CaRpgVthc 



{Vs4,d — V4,,g) 



(7) 



(8) 



' {Vatpd - Vk/2) ' 

Here Vacttd is the dust azimuthal velocity, vk is the local Ke- 
pler velocity, Cr is the dust-grain surface-area-to-mass ra- 
tio, pg is the local gas density (set equal to the value at 
the midplane in a vertically isothermal disc), and Wthcrm 
is the gas thermal velocity. We assume spherical grains of 
internal density 1 g cm~'^, and vary the grain radius, Sd. 
We consider grain sizes of 2sd = 0.2 fim - 2 cm, which is 
meant to represent the general range of grain sizes produced 
by coagulative grain-growth processes. Above some critical 
size, laboratory experiments and coagulation models show 
that bouncing and coll isional fragmentation represent barri- 
ers to further gro wth (JBrauer. DuUemond fc HenningJuOOa : 
IZsom et al.l 120101 ). The precise value of the critical size de- 
pends upon the gas disc conditions, and is not precisely 
known. We assume 2 cm as the upper-limit grain size for 
tracking a global dust distribution. Except in the case of 
instantaneous grain growth (discussed in ii l4.2|l . our simula- 
tion results are qualitatively identical for maximum assumed 
grain sizes as small as 0.2 mm. 

The local velocity for turbulent diffusion is. 



Vt 



urb = \/2Dp/At , 



(9) 

where Dj, is the local diffusiv i ty of t he particle ensemble. We 
use the lYoudin fc Lithwickl ((20071) scaling for the particle 
radial diffusivity. 



1 + ^T^ 
'(1 + ^1) 



(10) 



where Dg — vjSc is the gas-particle diffusivity, and 
Ts — S^Jk/ (CiiPgUthorm) is the normalized gas-drag stop- 
ping time. In general we assume a Schmidt number of 
Sc = V jDg = 1, but also run simulations in the flducial-disc 
model using Sc = 1/2 and 2. Equation (1101) has be en numeri- 
cally validated bv lCarballido. Bai fc Cuzzil (|201ll ') . However, 
note that for Ts ii 1, it deviates from the earlier common ex - 
pression given bv lCuzzi. Dobrovolskis fc ChampnevI (|l993h : 
Dp,CDC = Dg/ (1 + St), where St = Ts/tc, and Te is a nor- 
malized eddy turnover timescale. We have run a compari- 
son simulation of the evolving dust distribution within the 
fiducial disc model using this other dust-diffusivity param- 
eterization (with commonly used t^ = 1) and have found 
that the results are qualitatively identical to those using 
the parameterization in Equation (|lUp . As discussed in iJ31 
headwind-drag generally confines grains to regions of the 
disc with Ts < 1. 



2.2.2 Initial Particle Distribution and Dust Mass 
Allocation 

We assume that, at i = 0, the entire mass of disc solids is in 
small particles (0.2 ^m) whose radial distribution matches 
that of the disc gas. Since coagulation occurs rapidly for such 
small-size particles, to a first approximation we expect that 
the grain-size distribution will rapidly approach an equi- 
librium dictated by collision/fragmen tation processes (e.g., 
iBirnstiel. Ormel fc DuUemondl (|201ll )). In more detail, how- 
ever, the outer regions of a disc will experience a time-delay 
in the production of large grains relative to the inner, denser 
regions. This should translate into a delay before which rapid 
inspiral of large grains can occur from the outer disc. 

The grain trajectories calculated in our simulations are 
non-coagulative, considering a single grain-size for the full 
lifetime of a given run. In order to account for the delay 
in the appearance of larger grain sizes, we employ a pair of 
simplified grain-growth models (described in Appendix [Bj 
to define a radially dependent timescale for the appearance 
of larger grains. Larger grains are then initialized at time 
to (R) > within ten radial zones for which grain-growth 
time scales (at characteristic radii) have been calculated. 
This initialization of all grains larger than 0.2 p,Tn diame- 
ter occurs only once per zone, explicitly neglecting the colli- 
sional or coagulative regeneration and resupply of any grain- 
size population that may be removed due to transport ef- 
fects. However, as discussed in ^ limiting the appearance 
of large grains (especially in the inner disc where they would 
be most likely to regenerate) has the primary effect of retain- 
ing that dust mass against loss inward onto the parent star. 
This simplification therefore results in the highest dust-to- 
gas ratio in the inner disc achievable with our disc-evolution 
and dust-transport models. 

Our simplified grain-growth models consider coagula- 
tion via differential settling (in the rapid-upper-limit case 
in the absence of turbulence) and random motions. We do 
not consider growth via differential-radial drift, which, e.g.. 
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Figure 1. Timing for the initiation of larger grains within our 
fiducial simulations for each of our radial zones. These correspond 
to the more rapid of the two growth timescales produced by our 
cartoon grain-growth models, discussed in Appendix |B] Points 
not plotted beyond a given radial position indicate that grains 
are not grown to that size within the disc lifetime. 



IStepinski fc Valageaa (|l99q ). found to produce more rapid 
growth time-scales, nor directly via turbulent velocity dis- 
persions. Our assumed dust distributions, therefore, may 
represent a lower-limit on the rapidity and radial extent of 
large-grain initialization. In H4.2I we consider the upper limit 
cases where large-grain formation is a) instantaneous, and 
b) independent of disc location (not confined to the inner 
regions of the disc) in order to define the boundaries of the 
behavior of our evolving dust/gas distributions within the 
grain-growth-model parameter space. 

Figure[T]plots the grain-growth timing employed within 
our fiducial disc simulations. We run our two cartoon grain- 
growth models independently and initiate larger grain sizes 
according the most rapid growth timescale indicated. Due 
to disc evolution and depletion, sometimes only one (or nei- 
ther) model will produce grains of a given size within a given 
zone. Our ten zones have characteristic radii spaced logarith- 
mically from 0.5 to 1000 AU, and zone boundaries placed 
logarithmically between. We simulate transport of 0.2 nm- 
2cm-sized grains, spaced per decade in grain size, running 
an ensemble of 2 x 10* grains for each grain-size transport 
simulation. Grains initialized within a zone are placed ran- 
domly following the f = to gas-mass distribution within that 
zone. 

Each simulation particle tracks an assigned gas- 
equivalent mass, 771GE, of solid material within the disc. The 
gas-equivalent mass can be used to directly measure the lo- 
cal change in the dust-to-gas ratio (the dust enhancement 
factor, E) within the disc. By adding up all gas-equivalent 
mass within a radial bin, Mge ~ X^'^ge, and comparing 
that to the disc-gas mass within that bin, Mg, we can cal- 
culate the locally-measured dust-to-gas enhancement factor 
within our simulations. 

Mge 



E 



Mge 



M„ 



(11) 
1 indicates that the local disc metallicity 



is unchanged from its initial (solar) value, whereas E > 1 
indicates an enhancement in the local dust-to-gas ratio. 

Our simulations do not explicitly track refractory and 
icy material separately or include evaporation/condensation 
effects for grains crossing the snow-line. Instead, the com- 
position (and therefore relative mass of grains) is assumed 
to be that appropriate for the local disc-temperature condi- 
tions. While the enhancement factor is, therefore, the fun- 
damental output of our dust-transport simulations, it can 
be translated into a local value for the dust-grain surface 
density, Ep, separate from the disc gas. To do this, we solve 
for Sp at a given binned {R,t) — >■ {i,j) using those assumed 
compositions and the more general definition of the enhance- 
ment factor, 



E..-^ /^ 



(12) 



where {'Ep/'Eg)^ = Zq refers to the initial (solar) dust-gas 
composition of the disc, the basel i ne me tallicity. Zq used in 
this work is taken from iLoddersI (|2003l ) and is broken into 
two components: 



Zq = ZoZrcl (T) 



(13) 



where Zq = 0.0149 is the total fraction of condensable ma- 
terial thought to be present in a solar-composition gas, and 
Zici (T) is the fraction of that material believed to actually 
be condensed based on the local gas temperature. The values 
of Z^ci used are 



Z,^i = 0.3289 for 

Zroi = 0.7129 for 

Zrci ~ 1 for 



r > 182K , 

182K > T > 41K , 

T<41K, 



with transitions corresponding to the condensation of water 
ice (the snow line) and to the condensation of methane (and 
ammonia) ice. 

Because initially all disc solids are contained within 
0.2 ^m grains, at t = all 0.2 fim grains are assigned 
TTiGE = Mofi/'i X 10*. When a new grain size, 2s d^ is initial- 
ized within a radial zone, all ttige values for grains 0.2 /im 
- 2sd presently within that zone are reassigned to follow an 
updated mass allocation for the grain-size distribution. For 
this mass re-allocation, we assume a power-law distribution 
of grain sizes, 

£-^- (^^) 

where n^ is the number density of grains of a given radius 
Sd, and 5s is the power-law index of the distribution. For 
the fiducial simulations, we take qs — —3.5, correspond- 
ing to the size distribution of grain s measured for the ISM 
(jMathis. Rumpl fc Nordsiecklll977l ). However, observations 
linked to grain -growth within di scs often predict flatter val- 
ues of qs (e.g. iRicci et al.ll2011^ . and in |4] we also report 
results that include qs = —2.5 and -4, the later commonly 
used to dep i ct a p urely coUisional size distribution, e.g., 
iBaifc Ston^(|2010al l. 

To calculate the fraction of the total dust mass within a 
zone to be rea llocated into a given grain-size bin, we follow 
ICaraudl (|2007h . where, for a general value of g^, 
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Pv 



*+1s 



J+1 



111/ 



l*max *miii 



(15) 



where pj, is the local volume density of all solids, Pp^s^ is the 
volume density of solids within size bin Sd, Smin and Smax 
are the lower and upper bounds, respectively, on the entire 
size distribution, and Sbmin and Sbmax are the bounds on a 
given size bin. For qs — —4, 



Pp 



In (Sbn 



■ In (Sbmin)] 



In ( Smin)] 



(16) 



We define size bins logarithmically, so that, e.g., redistribut- 
ing mass across 0.2 /xm and 2 fim size bins {sd ~ 0.1 
and 1 ^m) has bin boundaries at Sd = 0.0316, 0.316, and 
3.16 ^m. Once we know the fraction of gas-equivalent mass 
to be portioned into each size bin, we divide it evenly among 
the particles of that size currently within the zone being up- 
dated. 

Note that because our dust-transport trajectories are 
non-coagulative, tracking only a single grain size per model 
run, grain growth in our simulations in each radial zone oc- 
curs only when our grain-growth models indicate that the 
next largest grain size should be initiated. At that time, all of 
the dust mass currently within that zone is summed and re- 
distributed according to Equations (|15|l or (|16|) , thereby up- 
dating the moE values of all particles (including those newly 
initiated) currently within that zone. We therefore main- 
tain a roughly equilibriated grain-size distribution within 
our disc during the early period of growth up to large-grain 
sizes. However, as radial drift proceeds and preferentially 
depletes mass from the larger grain-size bins, this equilibra- 
tion is not maintained. Furthermore, we assume that larger 
grains will be built from grains only slightly smaller than 
themselves. Therefore, in terms of mass reallocation, we do 
not initiate new, larger grains in a zone if radial drift has 
removed all grains of the next smallest size from that zone. 

Without a mechanism to replenish the large-dust popu- 
lation, our simulations therefore produce a conservative es- 
timate of the degree of dust-mass infall experienced within 
an evolving disk. As (discussed in ij4]) grain growth and the 
infall of large grains primarily leads to the draining of dust- 
mass from the disc, we expect that our simulations repre- 
sent rough upper-limits on local dust-to-gas enhancements 
in the inner/main disc at times greater than the earliest disc- 
evolution timescales (t > 0.5 Myr.) For possible alternatives 
at earlier times, see M.2\ 

Finally, we report the results for our dust-distribution 
simulations as a map of dust-to-gas enhancement, E, within 
radial and temporal space, {R,t). While our simulations 
report mass values for each radial grid-point (600 total) 
of our model disc and at each output time-step, usually 
5t — 200 years, we present our results at lower {R,t) res- 
olution in order to reduce noise from our particle-ensemble 
representation of the dust distribution. We present results 
using mass radially summed across three grid cells for 200 
radial grid-points of results output, and averaged across usu- 
ally 250 time-steps for a results At — 5 x lO'' years. 
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Figure 2. Comparison of dust-transport and grain-growth 
timescales for 0.2 mm grains within the fiducial disc model. 
Timescales for inward advection, inward drift, and difi^usion of 
grains are plotted for t = 10^ years of the fiducial disc model. 
Timescales for grain-growth are from t = and include disc evo- 
lution. There is a break in the plotted advection timescale where 
the dust flow is radially outward. The dashed line references the 
disc lifetime. 



3 RELEVANT TIMESCALES AND 
DISC-MODEL CONSTRAINTS 

We begin by considering several timescales relevant to grain 
transport, which are plotted in Figure [2] for 0.2 mm grains 
at early times within the fiducial disc model. The first is 
the time-scale to grow grains to a given size, for which we 
plot the initiation timing taken from the models shown in 
Figure [l] Next is the time-scale for grains to spiral inward 
onto the parent star. This can be approximated as the dis- 
tance to the star divided by the local inward-drift/drag ve- 
locity. However, there are actually two inward velocities to 
consider. The first is the head-wind drag velocity by itself, 
which is commonly used in such loss estimates applied to 
static disc models. The second is the drag-advection veloc- 
ity on the grain, calculated as Vsrd from Equations (O & (O 
above, which includes the (mostly) inward gas flow resul- 
tant from disc evolution. Within our model, the drift-only 
velocity, Ur,drift, may be calculated with the same method 
as Vsrd only setting Vr,g — 0. The two inward-loss timescales 
are therefore defined as 



idrift = —R/Vr, drift ■ 
f'advcct ^1 '^ srd • 



(17) 
(18) 



The final timescale to consider is the timescale for turbu- 
lent/diffusive redistribution of the particle ensemble, which 
for Figure [2] we define as 



= R'/Dj, . 



(19) 



From the timescales plotted in Figure (2] one sees that 
there are several important trends to consider. 

(i) The shortest time-scale is that for the appearance of 
these 0.2 mm grains within the disc, out to at least a few 
tens of AU. At any one radius, then, it would be reasonable 
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to assume that the initial phases of grain growth are effec- 
tively instantaneous. However, when we consider the disc as 
a whole, we find that even our lower limit on the growth 
time scale in the outer disc exceeds other time scales in the 
inner disc. This implies that by the time sub-mm grains 
have formed at tens of AU, we expect the major fraction of 
similar grains to have been lost from the inner few AU. 

(ii) There is a large difference between the time-scales 
tdrift and iadvect- The inward drift time scale, which would 
represent the flow of particles in an inviscid disc, is rela- 
tively long and almost constant with R. These features are 
consistent with the idea that radial drift results in a pile-up 
of solids in the inner disc. In an evolving disc, however, the 
advection time scale is substantially shorter everywhere in- 
side 10 AU. Solids can therefore be lost via advection along 
with the gas, without having had time to pile-up. 

(iii) Finally, the advection and diffusion timescales are 
similar, except in the outer disc where headwind-drag infall 
dominates. Turbulent redistribution of grains can therefore 
matter in the inner disc. Since turbulence represents a dif- 
fusive process, it will act to limit sharp radial gradients in 
particle concentration. 

For planetesimal formation, two basic properties of the 
background disc model are particularly important: the solid- 
gas coupling strength for particles of a given size, and the 
magnitude of the deviation from Keplerian rotation due to 
gas pressure gradients. These quantities enter into models 
for planetesimal formation via either streaming instabilities, 
or vertical precipitation followed by collapse. It is therefore 
useful to discuss how they vary in the gas disc model that 
we assume. 

Local dust-gas coupling is parameterized using the lo- 
cal gas-drag stopping time on particles within the disc. This 
depends on the size of the particle considered, but for char- 
acterizing the disc, we can plot the local particle-size corre- 
sponding to Ts = f^Kistop = 1, which is a stopping time of 
particular interest to some theories of planetesimal forma- 
tion. In the Epstein-drag regime 



6x10' 



0.2 mm 






(20) 



where Sd and pd are the dust-grain radius and internal den- 
sity, Pg is the local gas density, and iithorm is the thermal 
velocity of the gas particles. At the midplane in a vertically- 
isothermal disc model. Equation (I20|) translates to 



2E„ 



(21) 



In Figure [S] we plot the Ts = 1 contours for pd = 
1 g cm^'^ grains within the fiducial disc model, overlain 
on the 0.2 mm component of our simulated dust-grain dis- 
tribution. We have highlighted the contour corresponding 
to Ts = 1 for 0.2 mm grains, and one can see that the 
dust distribution within the disc is confined well interior 
{ts < 0.1) to this contour. This is consistent with the time- 
scales for radial drift presented in Figure [2] (ts = 1 occurs 
near 200 AU at the peak in the head- wind-drag effect), and 
also w ith the simulations of iBrauer. DuUemond fc Henninej 
((2003) , though radial confinement in our fiducial simulations 
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Figure 3. r^ = 1 grain-size contours {2sd = 0.2^in-2 m) overly- 
ing our simulated distribution of 0.2 mm grains within the fidu- 
cial disc model. Enhancement values of the distribution are scaled 
to the total simulated dust-grain distribution, of which 0.2 mm 
grains make up some fraction. Note that Ts = 1 contours follow 
Eg contours in a vertically isothermal disc (Eg. I21II . 



is about an order-of-magnitude greater than theirs, perhaps 
due in part to our approximate treatment of grain growth. 
We find identical behavior (confinement to Ta < 0.1 regions 
of the disc) for the simulated distributions of each dust-size 
component we have run, including only minor variation if 
the diffusivity (1/Sc) is varied by a factor of two. Headwind 
drag therefore confines larger dust solids to the inner/main 
disc. 

The pr ecipitation mechanism of planetesimal formation 
explored bv lYoudin fc Shul (|2002h is formally independent of 
the local dust size distribution (though grains do need to be 
large enough to settle toward the disc midplane, particularly 
if intrinsic disc turbulence is important for the vertical dust 
profile). However, streaming-instability models for planetes- 
imal formation find that particle clumping and collapse oc- 
curs specifically f or particles that are marginally coupled to 
the gas motions (IJohansen et al.l 120071 : Ijohansen fc YoudinI 
I2OO7I : iBai fc Stondl2010a^ . This means particles with nor- 
malized stopping ti mes, Ts, preferentially near 1 (though 
iBai fc Stond l|2010al ) report that grains as small as Ts ~ 10~^ 
still participate in the streaming instability). Because the 
headwind-drag effect dictates the rapid removal of grains 
from their Ts ~ 1 disc locations, it is hard to reconcile plan- 
etesimal formation models that require Ts ~ 1 with aero- 
dynamic grain-transport theory. Observationally, however, 
chondrules (~m m in size) do make up a large frac tion of 
most meteorites (|Cuzzi. Davis fc Dobrovolskisll2003 ). Plan- 
etesimals in the inner region of the Solar Nebula may then 
have been largely built from mid-to-large sized grains (see. 
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Figure 4. 11,, contours (red) overlying Ts = 1 grain-size contours 
(black) for the fiducial-disc model, n,; contours are for values of 
0.05, 0.1, 0.25, 0.5, and 1. Ts = 1 grain-size contours are spaced 
in decades 0.2 fim-2 m. 



e.g., ICuzzi. Hogan fc Bottkg (|201Cl ) for a discussion), and 
planetesimal formation models using such grain sizes may 
be well justified. 

Finally, a requirement of planetesimal formation 
common to both precipitation and streaming-instability 
models (and consistent wit h the model presented in 
ICuzzi. Hogan fc Bottkj (|201CI )) is that the local radial pres- 
sure gradient within the disc is small. The radial pressure 
gradient can be parameterized via. 



Vs,t> = 



{VK - V4,,g) 



Vk 



(22) 



since it is the source of the offset between gas and Keple- 
rian o rbital velocities. The dependence of the lYoudin fc Shul 
l|2002l ) precipitation mechanism on low rjsti, is explicitly out- 
lined in ^ Some papers reporting simulations of streaming- 
instability clumping of particles instead use fir, = rji^v-KJCa, 



n„ 



{VK-V^^g) 



(23) 



and in Figure |4] we plot If,, contours in red overlying the 
Ts = 1 g rain-size contours for the fiducial-disc model. In 
iBai fc Stone (2010b| ), the authors report simulations of par- 
ticle clumping via the streaming-instability in which they 
vary 11,, between 0.025, 0.05, and 0.1. They find clumping 
requires Ts JJ 10~^ for most particles and super-solar metal- 
licity when 11,, = 0.05, and a monotonic increase in required 
metallicity for increasing 11,,. Comparing these results to the 
spatial distribution of 11,, values within our disc model, it is 
clear that even the purely gaseous properties of the disc work 
to disfavor planetesimal formation at large radii. It remains 
to be explored whether potentially significant dust-to-gas 
enhancements above solar can be produced closer in. 



4 RESULTS: THE GLOBAL DISTRIBUTION 
OF SOLIDS 

4.1 Fiducial Results 

Figure [S] shows the simulated solids enhancement (change 
in the dust-to-gas ratio) as a function of R and t within 
the model disc. It is overlain with contours showing the 
evolution of the gas surface density. The bottom panels of 
the figure plot contributions from the individual grain-size 
components that make up the total solid mass. (Note that 
the green, scallop-shaped enhancement contours in the in- 
ner/main disc in the composite map are artifacts of the size- 
binning, matching the terminal-infall boundaries of succes- 
sively larger grain sizes.) In these maps, green areas repre- 
sent no change from the t = dust-to-gas ratio, red areas 
represent depletion of dust relative to the gas, and black 
areeis represent x2 enhancement of the dust or greater. 

From the figure, we see that the fiducial-disc simula- 
tions do not produce large radial concentrations of dust in 
the main and inner discs until near the end of the disc life- 
time. In ij5l we will compare these measured enhancement 
values to predicted requirements for the precipitation model 
of planetesimal formation. Recall, however, that a common 
minimum requirement quoted to aid planetesimal formation 
is a doubling of the dust-to-gas ratio above solar metallic- 
ity, and that, as discussed in ijl] it is the first few million 
years of disc evolution that are of greatest interest in terms 
of bulk planetesimal formation. Potentially significant en- 
hancements are not produced by our models during this 
phase of disc evolution. 

Figure [5] illustrates several other characteristics of the 
simulations. At early times, a noticeable portion of the 
small-grain population is advected outward as the disc ex- 
pands. However, inward drift of particles quickly sets in, and 
leads to the depletion of dust relative to the gas in the outer 
regions of the disc. In the inner disc, this influx of mass is 
in turn offset by concurrent loss (due to further inspiral) 
of most of the larger-grain population. The net result is a 
significant loss of dust mass from the evolving disc. 

The significance of the mass contained within the larger 
grain sizes is depicted in Figure |6] where, for t = 1.5 Myr, we 
plot the enhancement factor produced when the results of 
our fiducial simulations are compiled (using sub-sets of the 
simulation output) for grain-size distributions truncated at 
successively smaller particle sizes. Here, we see that grain- 
growth to roughly 20 /im sizes results in inward drift that 
depletes the outer disc of dust, and it is for grain-size trun- 
cation at 20/im that our simulations produce peak enhance- 
ment ratios within the inner/main disc. Growth to larger 
sizes, however, subsequently depletes the inner disc of this 
dust-mass excess, bringing enhancement values back toward 
unity. 

Finally, we can solve for a representation of the local 
surface-density of dust particles, Ep, using Equations (|12[) 
and (|13p in order to gain a sense of the diminution of solids 
available in an evolving, accreting disc as time progresses. In 
Figure [T] we plot the time evolution of Ep and of the local 
dust-mass per AU within the fiducial-disc simulations. The 
snow-line is clearly marked by jumps in the distributions at 
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Figure 5. Composite map of dust/gas enhancement factor (top panel) within the fiducial disc model and maps of some of the contributing 
grain-size bins (bottom panels). Regions colored in black have a solids/gas enhancement factor of x2 or greater. The black contours trace 
the evolving gas surface density, Eg, spaced by orders of magnitude between 10^'* and 10'^ g cm~^. The grey regions mask out bins for 
which the particle statistics were too low to report dust/gas enhancement values. 
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Figure 6. Measured enhancement of local dust/gas for the 
fiducial-disc model simulations at t = 1.5 Myr varying the upper- 
limit on the grain-size distribution. For no grain growth (0.2 fim 
grains only) the dust enhancement factor remains everywhere 
near unity. Allowing larger grain sizes allows inward drift and 
some pile-up of the dust population. 



around a few AU. At i = 0, our fiducial disc lias just over 
100 Earth-masses of solid material spread throughout the 
entire disc, with the greatest concentration just beyond the 
snow line at around 5-10 AU. As the disc evolves, the solids 
distribution spreads outward with the disc and decreases in 
total quantity, to 40 i\fEarth at t = 5 x 10^ years, and 10 
MEarth at t = 2.5 Myr. 



4.2 Varying Disc Model Parameters 

To quantify the role of radial dependence on dust/gas en- 
hancement results through both the local maximum of grain 
sizes considered and the time delay in the formation of the 
largest grains, we present two simulations using alterna- 
tive grain-growth model assumptions. In the first, the initial 
particle-size distribution does not vary as a function of disc 
radius. All particles of all sizes (0.2 fira ~ 2 cm) are instead 
initiated at i = evenly distributed throughout the entire 
disc-gas-mass distribution. This is similar to the radially in- 
de pendent size distribut ion used as an initial condition in 
the lYoudin fc ShrJ (|2002l ) calculations of dust pile-up via ra- 
dial drift, except that we retain the simple power- law form 
of the size distribution with Qs = —3.5 throughout, while 
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Figure 8. Enhancement map for a simulation using a radially 
independent initial grain-size distribution. 



Figure 7. Time evolution of (a) the dust surface-density profile, 
and (b) the dust mass-per-radius distributions for the fiducial-disc 
simulations. 



lYoudin fc ShtJ ((2003) considered a two-component distribu- 
tion representing chondrules and matrix material. 

The results of our simulations using a radially indepen- 
dent grain-size distribution are plotted in Figure [8] and are 
distinct from those of our fiducial simulations. Due to the 
large fraction of total dust mass now residing in large grain 
sizes, most of the dust-mass in the disc rapidly falls inward 
to small AU. It is thus quickly lost onto the parent star, 
but not before creating a period of distinct enhancement 
in the dust-to-gas ratio in the inner disc. The peak in this 
enhancement occurs at i ~ 4 x 10'' years at about x5 so- 
lar (initial) metallicity. (If we instead assume a maximum 
grain size of only 2 mm, the peak in the enhancement oc- 
curs at i ~ 2 X 10^ years.) After this peak, disc metallicity 
declines swiftly, and by f = 5 x 10''' years the metallicity 
across the whole disc is at a fraction its f = value. While 
this extreme-case grain-size distribution does produce sig- 
nificant enhancements in the dust-to-gas ratio early in the 
disc lifetime, the time to reach this enhancement peak is 
shorter even than the expected time-scale for disc assem- 
bly (~ lO' years). It is thus hard to envisage a disc forma- 
tion scenario in which the pile-up seen in this instant-grain- 
growth-everywhere simulation could be physically realized. 
Moreover, such a short window supplied by this simulation 
for planetesimal formation is in direct contradiction with 
meteoritical evidence implying large-body formation in the 
Solar System persisted for millions of years. 



In the second alternative simulation, we again assume 
that grain growth is instantaneous, initiating all particles at 
i = 0. However, we now return to the assumption that the 
largest grains will grow only in the inner, dense regions of 
the disc. The results depicted in Figure [9] use the same sim- 
ulation output as those of Figure |8] but have been compiled 
selecting only those grains initiated in the zones appropriate 
for their size as defined by our cartoon grain-growth model, 
e.g., 2 mm grains initiated at i? < 9.6 AU.[jThis simulation 
produces results similar to those of our fiducial simulations 
but with two notable exceptions. First, the small-grain pop- 
ulation is depleted relative to the fiducial simulations, pro- 
ducing a somewhat greater depletion of dust relative to the 
gas at large distances and at late times. In the fiducial sim- 
ulations, small grains have time to advect and diffuse into 
the massive outer disc before dust mass at smaller distances 
is converted into larger sized grains. 

Second, there is a peak in the enhancement factor of 
just over x2 solar metallicity at small radii around t = 
6 X 10'' years. This peak corresponds to the terminal infall of 
the 0.2 mm sized grains, which are carrying a larger portion 
of the dust mass than in the fiducial simulations. Assuming 
instantaneous grain growth also allows infall of these larger 
grains without delay from the further regions of the disc 
(out to i? = 52 AU). While the 2 mm and 2 cm grain-size 
bins are also holding more dust mass than in the fiducial 
simulation and infall more swiftly than the 0.2 mm grains, 
their presence is too radially confined to produce an earlier 
peak in the dust/gas enhancement like that seen in Figure [H] 

^ Note: To compare to Figure[T] 9.6 AU is the outer boundary of 
the radial zone with characteristic radius defined at 6.3 AU. 



© 0000 RAS, MNRAS 000, 000-000 



Particle evolution in protoplanetary discs 11 



6x10 




6x10 



Low Diffusivity (Sc=2) 



0.0 0.5 1.0 1.5 2.0 
Enhancement factor 



Figure 9. Enhancement map for a simulation witti all grains 
initiated at t = 0, but with maximum grain sizes set by radial 
zone using our cartoon grain-growth model. 



The enhancement peak produced in this simulation is timed 
more favorably to assist in planetesimal formation than the 
fully radially independent simulation, particularly under the 
potentially lower enhancement requirements of streaming- 
instability-type planetesimal formation. However, its dura- 
tion is still noticeably shorter than the meteoritical evidence 
would seem to indicate necessary. 

The results of these simulations therefore seem to indi- 
cate that faster grain growth leads to a greater probability 
of dust pile-up in the inner regions of the disc, but that 
this effect is limited by the mass supplied in large grains 
initially a large distances from the star. Furthermore, faster 
growth may still be insufficient to supply an extended period 
of dust/gas enhancement. For the remainder of our simula- 
tions, we return to initiating particles radially and tempo- 
rally based on our cartoon models of grain growth. 

Figure[TO]plots the enhancement distributions produced 
varying the Schmidt number, Sc = v/Dg (inversely varying 
the gas-disc scaling for the particle diffusivity), by a fac- 
tor of two in either direction. We find the precise value of 
Sc used to be unimportant to our results. A lower diffu- 
sivity leads to greater segregation of the particle-size distri- 
bution and emphasizes the discrete particle sizes used for 
these simulations. However, a real disc should have a mostly 
continuous particle-size distribution, and therefore the re- 
sults for the global distribution of solids cannot be said 
to be qualitatively different from the fiducial model. The 
higher-diffusivity case allows large grains to remain in the 
disc slightly longer than usual (all 20 /im grains are lost 
after 4.1 Myr rather than by 3.8 Myr as in the fiducial sim- 
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Figure 10. Enhancement maps from simulations run with Sc = 
1/2 and 2. We use the fluid-dynamics convention that Sc = u/Dg, 
while retaining the local scaling for particle difi^usivity based on 
grain size using Equation jTo}. These simulations were run with 
output at lower resolution than in the fiducial setup and therefore 
compiled using six radial-grid-points per radial-bin and averaging 
25 time-points for every At = 5 X 10* years time-bin. 



ulations) , but again, this is insufficient to substantially alter 
the global distribution of solids. 

In Figure [TT] we plot the composite enhancement maps 
for simulations run varying three disc-model parameters: the 
initial disc mass, Md,o', the initial disc compactness, Rd, 
and the magnitude of the disc viscosity, a. Qualitatively, 
the patterns of evolving global enhancement are roughly 
constant across the different disc-model setups. The infall 
of the smaller grain sizes causes depletion in the dust rel- 
ative to the gas at the outer edge of the disc, while in the 
main disc enhancement values remain near unity for some- 
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Figure 11. Enhancement maps for simulations run in several different disc models, separately varying Rd (horizontal), Mnfi (vertical), 
and a (diagonal). Fiducial-model map {M^) g = O.O3M0, R^ = 20 AU, a. = 10~^) showrn in the center panel. Note that vertical scales 
are each set to the same years-per-inch to illustrate the varying rates of disk evolution. 



where between 40 and 60% of the total disc lifetime. The 
initially more compact disc shows somewhat greater outer- 
disc depletion because a larger fraction of grains are pro- 
cessed to larger sizes at early times, while the disc using 
a = 5 X 10""^ shows some of the enhanced segregation of our 
discrete particle sizes typical of the lower-diffusivity sim- 
ulation in Figure [lO] (as, in our transport model, lower a 
means lower disc viscosity and diffusivity). However, the ba- 
sic shape of dust-to-gas values within the disc as controlled 
by aerodynamic forces appears fairly universal. The primary 
result, of these simulations then, is that a smoothly-defined, 
azimuthally-symmetric, evolving disc model tends to lose 



dust-mass smoothly inward onto the parent star at a rate 
that generally keeps pace with the loss of accreting disc-gas 
mass, and that variety in basic disc properties/parameters 
does not tend to produce special-case discs with respect to 
this process. 



5 COMPARISON WITH REQUIRED 

ENHANCEMENT FACTORS 

In this section, we compare o ur simulated enhancem ent fac- 
tors to those required by the lYoudin fc Shul (|2002l ) precipi- 
tation model of planetesimal formation. In the precipitation 
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model, planetesimal formation is accomplished by the grav- 
itational collap se of so li ds set tled out to the disc midplane. 
As outlined in ISekival (|l99a ), it assumes that the disc is 
quiescent, with no global turbulence to stir the particles up- 
ward, and that grains will settle out to the limit provided by 
the Kelvin-Helmholz shearing instability. For midplane dust 
densities in this settled state above a certain mass threshold, 
the gravity of the dust sub-disc should dominate and lead to 
the collapse of so lids at the midplane. The criteria for this 
collapse, given in lYoudin fc Shul (|2002l ) Equation (15), is 



/RicVs<pRpgS {ip) 



(24) 



where Ep_c is the critical surface density of the solids for 
collapse, Ric is the critical Richardson number defining the 
degree of dust settling prior to the onset of the Kelvin- 
Helmholtz instability, and 5(1/)) is a correction factor ac - 
counting for the self-gravity of the gas. lYoudin fc Shul (|2002l ) 
define 



{dP/dR) 



'Id 4) — 


2pgRQ? ' 


and 




sW 


= (1 + V)ln 


^ 


-ilvGpg 



(25) 



{l + i> + ^l + 2'4>) 



\/i + 2V; 



(26) 



For a gas density of pg ~ T,g/2Hg, where Hg — c^ R/vk in 
a vertically isothermal disc model: 



(Ep/E 



a)Q 



- VK Vs,ps (V;) 

' Cs ZoZrcl (T) 



(27) 



using the conventions for Zo and Zici (T) defined in ii2.2.2l 
One expects that any degree of turbulence beyond the 
Kelvin-Helmholz effect should stir the particle distribution 
to greater vertical heights than those assumed in this Ric 
scenario; certainly, the a = W~^ used in our fiducial disc- 
evolution models suggests that a significant fraction of the 
small dust grains should be vertically well-mixed with the 
gas. However, assuming that grains could settle out to their 
Kelvin-Helmholtz limit allows us to calculate the absolute 
minimum enhancement-factor required for precipitation of 
planetesimals to proceed in that assumed non-turbulent en- 
vironment. 

The precipitation-collapse criterion is separate from the 
simpler gravitational-coll apse criterion u s ing th e Toomre-Q 
parameter, Qp < 1. From lYoudin fc Shul (|2002l '). 



H„ 



TrGrEp.c 

VRicVStpRhd') 



(28) 



h{^p) = \/l + 2?/> - ?/; In 



(1 -n/> -I- yi -I- 2^) 



The minimum enhancement criterion for Q-coUapse is then 

VKVs^h (Tp) 



-7vGRJ^gZoZ,,i{T) 



(29) 



model Sau (g cm ^) -qs Tau (K) -qt 



H 


1700 


3/2 


280 


1/2 


A 


1700 


3/2 


170 


0.63 


B 


1700 


3/2 


100 


3/4 


Af 


1700 


1 


170 


0.63 


Af5 


8500 


1 


170 


0.63 


fiducial 


1180 


1 


280 


0.6^0.5 



Table 1. Parameters for disc models used in JYoudin &: Shul 
1I2OO2I ) versus rough equivalents for our fiducial disc near 
t = 0. In general, Eg (R) = Sau (R/IAU)''^ and T (R) = 
'Tau (^/lAU)''^ . "H" de signates the mi nimum-Solar-Nebula 
model of Hayashi-1981 (see JArmitagei | J2010|') pp 4-5). 



which tends to dominate in the outer regions of massive, 
cold discs. 

While the critical Richardson number is com monly set 
to Ri c = 1/4, recent numerical simulations by I Lee et al.l 
(J2OI0I ) find the onset of Kelvin-Helmholtz instability occur- 
ring at 



Ric' 



£p,0 

36 



(30) 



where tp.o is the current dust-to-gas volume-density ra- 
tio at the disc midplane. ep,o can be written as ep,o « 
(Ep/Eg) / (Hp/Hg) and, if we continue to pursue to limit of 
maximum settling so that Hp follows from Hp^c in Equation- 
set (|28p . we can write Ric from Equation (|30|l in terms of 
the current dust-to-gas enhancement and a Ri value defining 
the current degree of settling. Absolute maximum settling 
should therefore occur to ffir("Eci30)l= Ri{'E(^^= Ri'c, or 



Ri'c = 






2/3 



(31) 



Below, we compare our simulated dust enhancements to pre- 
cipitation/collapse requirements using Ric — Ri'c, so that 
Equations ^} & ([29]) become 



and 



Eq,> 



6c,^O^rel(T)0r(VO ' 



VKr)s<i,h (tp) f r^K 



6ZoZ,ci(T) \-kGE 



3/2 



(32) 



(33) 



This choice results in a slight decrease (up to a factor of 
two less) in minimum dust enhancements required in the 
inner/ main disc to produce planetesimal formation relative 
to the lYoudin fc Shul (|2002h calculations, which use Ric = 
1/4. 

Figure [12] plots the req uired enhancement fa ctors for 
the disc models discussed in lYoudin fc Shul (|2002h (panel- 
a) and for our fiducial disc at several times (panel-b). In 
Ta ble [1] we list the ma in parameters defining the differ- 
ent |Youdin_&_Sliuj (I2OO2I ) disc models, as well as the values 
roughly correspond ing to the inner regi ons of our fiducial 
disc at early times. (|Youdin fc Shul (|2002l ') use simple power- 
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Figure 12. Requi red Enhan c emen ts for midplane precipitation 
using the results of lLee et alj 1120101) and Eq uation II31I I to calcu- 
late Ric for the discs of lYoudin &: Shul 1 120021) and our fiducial-disc 
model. In our fiducial disc at t = 0, the snow-line minimum value 

IS -C/TircciD ^ o.oo. 



law disc models while our models have an exponential tail-off 
of surface-density in the outer disc.) 

Figure [12] demonstrates several interesting features of 
these enhancement requirements. 

• First, note that our fiducial model calls for remarkably 
high -Eprecip in the outer disc, particularly at early times. 
This is due to the strong dependence of -Eprocip on rjs^,. When 
the radial pressure gradient is steep, as near the outwardly 
expanding disc edge, shearing between the dust and gas is 
strong, requiring large dust-to-gas ratios to be overcome. 
This limiting requirement is shared by some other proposed 
models of planetesimal formation, such as the streaming- 
instability model, as disc ussed in ^ 

• Second, while the lYoudin fc Shul ()2002l 'l paper de- 
emphasizes the roles of the disc gravity, represented by s (tp), 
the time-series for -Eprocip of our fiducial model, particularly 



in the inner disc, demonstrates the importance of the local 
disc mass for meeting dust/gas-precipitation requirements. 
s (-(/)) (and therefore -Eprecip) becomes large, when tp is small, 
which may occur for either a low local-disc surface density, 
or for a locally hot disc. 

• Finally, the -Eprocip curves for both sets of disc models 
emphasize the importance of the ice lines, particularly the 
snow line, as places of potential local minimum in metal- 
licity requirements for planetesimal formation. While our 
particle-transport model does not treat e vaporation and 
vapor -transport effects, the simulations of ICiesla fc Cuzzil 
()2006l ) do show an additional peak in the dust population 
just outside the snow line due, at least partially, to local 
concentration of diffused and recondensed vapor. This also 
corresponds to a peak in the disc opacity and heating, how- 
ever, so it is not clear what the degree of advantage the 
region just outside the snowline lends to planetesimal for- 
mation via gravitational precipitation. 

In Figure 1131 we compare minimum enhancement fac- 
tors required for planetesimal precipitation to enhancement 
results from our simulations compiled using three different 
assumed particle-size distributions. From this figure, it is 
clear that the size distribution used to compile our sim- 
ulation results has little impact on the conclusions of the 
fiducial runs. Radial drift of particles does not produce en- 
hancement sufficient to lead to planetesimal formation via 
precipitation for this model disc. The exception occurs at the 
very end of the disc lifetime when photoevaporative clearing 
of the disc leads to an outward-sweeping pressure-maximum 
point at large AU where rjg^ goes to zero and -Eprccip drops 
to small values. 

The first two panels of Figure [13] also include the results 
of the simulation discussed in > j4.2l where all particle sizes 
were initiated at i = and distributed throughout the en- 
tire disc gas mass, with no radially-dependent grain-growth 
restrictions. While this simulation yields a distinct peak in 
the dust-to-gas enhancement early on in the inner disc, it 
is still less than the minimum enhancement required by the 
precipitation model. 

Finally, while the simulated enhancement distributions 
vary little between different evolving-disc models (as shown 
in Figure llip the enhancement factors required by the pre- 
cipitation model should vary according to factors such as 
disc mass. In Figure 1141 we plot the enhancement criteria 
for precipitation near the beginning and end of the disc 
lifetime for the various discs used in generating Figure 1111 
While more massive discs, in particular, produce lower mini- 
mum enhancement requirements, these requirements remain 
generally high relative to the dust-to-gas enhancements pro- 
duced in our simulations. 



6 COMPARISON TO DISC OBSERVATIONS 

One of the most general results from our simulated dust dis- 
tributions is that we predict the outer-most edges of a disc 
should be depleted in dust relative to the gas. Furthermore, 
we predict that the dust that does reside in the outer reaches 
of the disc should be fine-grained in order to remain mixed 
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Figure 13. Measured enhancement values for simulations in the fiducial-disc model at several times and using several particle-size 
distributions. Compared to the enhancement factor required for precipitation/collapse as given by Equations II32I I &: I I33I I. The dashed 
curves plot the simulated distribution using the to = 0, iJ-independent grain-size distribution discussed in i|4.2l which in the left-top-most 
panel is plotted at the peak in enhancement, t = 4 X 10* years. 



with such tenuous gas. It is not yet clear from observations 
whether the outermost regions of a disc tend to be gas rich, 
though our prediction is consistent with some observations 
that find smaller disc radii w hen examining the d ust versus 
the gas signatures of a disc (JHughes et al.l 120081 ). As good 
observations at the outer edges of a disc require high sensi- 
tivity, a full test of outer-edge disc metallicities awaits more 
sophisticated facilities, such as ALMA. 

However, there are observations of dust grain-size dis- 
tributions within discs around some stars that are incon- 
sistent with the dust-distribution simulations we have pre- 



sented. Our simulations predict a global loss of mm-sized 
and larger grains from the disc within half a million years 
(even in our simulation of H4.2l where grains of all sizes were 
placed out to the largest distances within the disc). How- 
ever, observations in the (sub-)mm bands see evidence for 
mm-sized grains, not just wit hin discs, but out to 100 AU 
or fart her within those discs (|Testi et al.l 120031 : iRicci et al.l 
l2010bl l . The cartoon grain-growth model we use to constrain 
the radial extent of simulation particles does not allow mm- 
sized grains to be initiated beyond 22 AU an any of the 
disc models we've run, and the inward flow from both accre- 
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tion and headwind drag certainly precludes any simulated 
mm grains fro m reaching 100 AU. Further more, even in the 
simulations of iHughes fc Armitagel ((20101) that considered 
grain transport (using the same fiducial disc model) within 
an outward flow of disc gas at the midplane, less than 6% of 
mm-sized grains initiated between 0.5 and 10 AU reached 
even 25 AU, and then only for a brief period of time before 
headwind drag forced them inward once more. 

That headwind drag tends to bar large grains from 
the outer dis c, in contradiction wi t h observations, i s well- 
known (e.g. IWeidenschillingI [l977l: iTakeuchi fc LinI l2005l : 



iBrauer et al.1 120071 ). iBrauer et al.1 ( 20071 ) specifically focus 
on attempting to solve the radial-drift problem by including 
increasingly sophisticated physics into their transport cal- 
culations, including dust settling and collective effects, and 
angular momentum exchange between a dust sub-layer and 
the main gas disc. The authors find these effects in general 
insufficient to slow radial drift of mm-sized grains to veloci- 
ties compatible with their observed persistence in discs that 
are generally 1 Myr old or older. 

The apparent inability of dust transport models to 
match observations poses a puzzle, since the basic one- 
dimensional description of the gas disc evolution appears to 
be a reasonably good fit to statistical measures of the popu- 
lation such as disc sizes, lifetimes, and accretion rates. It is 
possible that the axisymmetric gas disc model, while glob- 
ally descriptive, is qualitatively inadequate as a background 
model for studies of particle transport. One way in which 
this could happen is if local or non-axisymmetric structures 
are integral to accurately describing the large-scale distribu- 
tion and transport of large dust grains. Otherwise, explain- 
ing outer-disc, large-grain populations with current aerody- 
namic disc theory seems to require such extreme scenarios 
as a coUisional population of outer-disc planetesimals (of un- 
known origin), or else very massive outer gas discs of very 
low metallicity (due to the sustained depletion of large grains 



observed to be in residence) and currently unknown impact 
on global disc evolution. 



7 DISCUSSION AND CONCLUSIONS 

In this paper, we have presented simulations of the global 
(re)distribution of dust within protoplanetary discs. Our 
models include evolution of the gas-disc surface-density pro- 
file, aerodynamic advection and diffusion of the dust-particle 
ensemble within the gas, and simulated grain-growth con- 
straints confining the appearance of the largest grain sizes 
to the inner/main regions of the disc. We find that the global 
distribution of solids within a ID evolving disc model follows 
a fairly common pattern, and does not depend strongly on 
the values of parameters characterizing the initial disc, or 
on those describing the turbulence within it. Specifically, 

(i) Large, mm or cm-sized particles are lost rapidly onto 
the the parent star (within 0.5 Myr). 

(ii) Growth of dust grains up to tens-of-microns size and 
larger leads to the depletion of dust relative to the gas to- 
ward the outer regions of the disc (beyond ~100 AU). 

(iii) The dust-to-gas ratio within the main and inner discs 
remains near solar (its initial value) for at least the first 40- 
60% of the disc lifetime. 

In general, grains of a given size are confined to regions of the 
disc interior to at least Tb < 10~^ and experience accelerated 
infall as Ts (oc Eg) contours near the outer boundary of that 
dust-p opulation distri b ution evolve inward. 

In lYoudin fc Shul (|2002h , the authors present calcula- 
tions of dust-to-gas enhancement factors resulting from ra- 
dial drift of dust grains within a static disc, and find a 
steady increase in dust-to-gas ratios at small disc radii with 
time. This is in marked contrast to our own simulation re- 
sults, which do not support the radial-drift hypothesis as 
a means to increase the metallicity in the inner disc in 
time for bulk planetesimal formation. The difference be- 
tween the two calculations is primarily the difference be- 
tween a purely static- and an evolving-disc scenario. As out- 
lined by the drift/advection timescales presented in Figure[2l 
a static disc produces ever-slower inward drift as grains ap- 
proach the central star, providing a natural environment for 
grain pileup. An evolving disc, however, will advect mate- 
rial steadily inward onto the parent star. Furthermore, as 
the disc evolves and thins, drift due to headwind drag is 
accelerated, thus limiting the slowing of inward drift possi- 
ble, even presuming a zone of the disc were to be free of an 
accretion/advection flow. 

While our simulations do not paint a picture of large- 
scale increases to disc metallicity within an evolving proto- 
planetary disc, they do suggest a state of long-term continu- 
ity of conditions, perhaps for the first several million years of 
disc lifetime. This is in keeping with observatio ns that sug- 
gest the processing of sol ids in discs starts early (|Kwon et al.1 
I2OO9I : iRicci et al.1 l2010al ) but extends late, as evidenced by 
the spread i n ages measured for chondrules embedded within 
meteorites (jAmelin et al.ll2002l ). Furthermore, that our mod- 
els predict similar conditions across a range in disc-model 
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parameter space is in keeping with the apparent ubiquity 
of large-body processing within discs around other stars. 
The fraction of stars observed to have circumstellar debris 
discs (with scattered-light levels around a thousand times 
brighter than that of our own zodiacal dust and Kuiper 
belt) peaks at around 50% for 20 Myr old B and A stars 
ICurrie. Plavchan fc KenvonlbOO^ I. 

We have failed to find circumstances in which global re- 
distribution of solids, due to radial drift, results in substan- 
tial enhancements in the dust-to-gas ratio anywhere within 
an evolving disc. If large enhancements are in fact a pre- 
requisite for planetesimal formation, it may instead be nec- 
essary to appeal to local or non-axisymmetric structures to 
provide them, although the effect of such structures may 
also be constrained by the headwind-drag produced lack of 
Ts ~ 1 particles, discussed in iJS] Local disc structures that 
could prove important for dust concentration and/or trans- 
port include disc-opacity transitions, such as at the snow- 
line, spiral arms prod uced within gravitationally unstable 
discs (|Rice et ahluOOa ). and vortices within the flow of disc 
gas. In discs subject to the magnetorotational instability, 
some simulations suggest that MRI turbulence produces lo- 
cal s tructure sufficiently long lived to co ncentrate dust parti- 
cles (jJohansen. Youdin fc Klahii2009al ). Determining which, 
if any, of these processes is important for particle transport 
and planetesimal formation will require an understanding 
not just of the strength of disc turbulence, but also of its 
detailed physical nature. 
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APPENDIX A: CALCULATING THE 
EVOLVING DISC-TEMPERATURE 
DISTRIBUTION 

The model-disc temperature is calculated using a power-law 
parameterization of the energy-balanced temperature at the 
disc midplane. To do this, the disc is first evolved once using 
an iterative solver to calculate the energy-balanced temper- 
ature at each grid point and each time-step due to passive 
and viscous heating and radiative cooling. We then make a 
power-law fit from the outputted temperature history. 

For the energy-balanced temperature, we calculate 
three components to the disc heating. The flux of energy 
due to viscous dissipation near the disc midplane is given 
by 



fvi 



RuT^a 



dR 



(Al) 



where the viscosity, v cc T, is taken from the alpha- 
parameterization described in i]2.1l The flux due to heating 
via stellar illumination is 

i^etar = 2^3 eT^ = ^^ , (A2) 

where (jb is the Stefan-Boltzman constant, Tirr is the passive 
temperature due to stellar radiation, and the radiative effi- 
ciency is assumed to be e = 1. For our disc evolution models 
we use Li, — 5I/soiar and set the incidence angle to the disc 
at = 0.05 everywhere. Finally, heating due to background 
radiation is given by 



-Fcloud — So-BeTt-iQmj 



(A3) 



where Tdoud = lOK is the background temperature of the 
cloud environment. 

These are used to solve for the midplane tempera- 
ture of the disc follow i ng th e approximation derived by 
iNakamoto & Nakagawal (| 19941 '): 



^sT^i, 



\R^^^9 



-dR [r''+2^ 



where Tmid is the temperature at the disc midplane, tr = 
Kit'Eg is the Rosseland mean optical depth, and rp is the 
Planck mean optical depth. For the op acity, kr, we use the 
piecewise parameterization provided bv lBell & LinI (| 19941 ) in 
which kr depends on both temperature and the mean gas 
density (for which we use the density at the disc midplane. 
Pa, mid). For simplicity we assume tp/tr = 1, but note that 
the specific ratio chosen has negligible impact on the energy- 
balanced temperature calculated. The iterative solver seeks 
a convergence in Tmid according to Equation (|A4|) . updat- 
ing the corresponding values for the midplane density (for a 
vertically isothermal disc profile), opacity, and disc viscosity 
at each pass. 

The primary focus of this paper is the non-uniform na- 
ture of the dust/gas distribution in an evolving disc, an effect 
which is disregarded for the purposes of the energy-balanced 
disc-temperature calculation. In this model of disc temper- 
ature, the effects of the dust distribution are felt only in the 
amount of energy deposited via disc accretion. If growth of 
the dust grains were taken into account, dust opacity should 
drop for larger grain sizes, thereby allowing for less disc heat- 
ing via accretion. However, at the time-scales / masses and 
accretion rates considered in our models, the evolution of 
disc surface-density contours is reasonably convergent using 
either a static (cooler, passive disc) or our evolving-energy- 
balanced temperature distribution. If, on the other hand, 
massive infall of the dust population led to significant dust 
enhancement in the inner disc, dust opacities in the inner 
disc should increase. However, at early times, such infall is 
comprised solely of large dust grains and so we expect that 
the effect should be negligible. Infall of the small-grain pop- 
ulation occurs only when the disc is too thin for accretional 
heating to be a significant contributor to the disc tempera- 
ture. 

Next, we fit the temperature history of a given disc 
model to a parameterization following the form 



T{R,t)^TAu{t) 



R 

lAU 



grit) 



(A5) 



First Tau and qt are calculated at each time-step using the 
intersection at R = 0.2 and 200 AU with the outputted, 
energy-balanced temperatures. Next, a quasi-exponential 
decay in time is fit separately to both Tau (i) and qr (t) 
following the form 



x{t) = {xo — a;oo)exp 



(A6) 



-I- (TB (Tirr + rcloud) 



(A4) 



with Xo and a:oo representative of the initial and final states 
of the temperature evolution, and t^ and bx fitting constants 
chosen roughly by inspection. The parameterization con- 
stants used for the temperature evolution of each of the disc 
models of this paper are given below in Table IXT] These pa- 
rameterizations are used to calculate disc temperature (and 
hence, the surface-density evolution of the disc gas) during 
each of the dust-transport runs (one per grain size per disc 
model) used to build our dust/gas-evolution simulations. 
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Md,o (Mq) Rd (AU) 



Tau,o (K) Tau,oo (K) 

<?T,0 iT.oo 



^Tau 



(yr) 



Oo'-r 



0.01 


20 


10-2 


458 


280 


1.9 X 10^ 


1.03 








-0.59 


-0.5 


2.2 X 10^ 


1.09 


0.03 


5 


10-2 


500 


280 


2.4 X 10^ 


0.966 








-0.61 


-0.5 


2.7 X 10^ 


1.01 


0.03 


20 


2 X 10-2 


500 


280 


2.3 X 10^ 


0.917 








-0.61 


-0.5 


2.6 X 10^ 


0.959 


0.03 


20 


10-2 


500 


280 


3.5x10^ 


0.925 








-0.61 


-0.5 


4.0x10^ 


0.966 


0.03 


20 


5 X 10-3 


500 


280 


5.1 X 10^ 


0.907 








-0.61 


-0.5 


5.8 X 10^ 


0.943 


0.03 


40 


10-2 


472 


280 


5.1 X 10^ 


1.01 








-0.6 


-0.5 


5.9 X 10^ 


1.07 


0.09 


20 


10-2 


600 


280 


7.3 X 10^ 


1.00 








-0.65 


-0.5 


8.1 X 10^ 


1.03 



Table Al. Disc parameters and power-law temperature fitting constants (Equations IIA5I I & IIA6I I ) for the disc models used in this 
chapter. The t = temperature is lower if there is initially much less disc mass at small AU. The time-scaling constants, tT^u ^^'^ ^it 
give a sense of the relative lifetimes of the various disc models, though the lifetimes are more than an order-of-magnitudc longer than 
either of the scaling-constants. The parameters of the fiducial model are in bold. 



APPENDIX B: TWO CARTOON 
GRAIN-GROWTH MODELS 

Within the particle-transport simulations, we use two 
simple-case grain-growth models to place constraints on 
when and where grains of a given size will appear within 
the model disc. These models consider growth of a particle 
due to differential settling from some height within the disc 
(called the raindrop model because it mirrors the growth 
of raindrops falling through the atmosphere), and particle 
growth due to random, thermal motions at the disc mid- 
plane. The basic schematic for these two models is given in 
Figure [BT] 

The raindrop mod el is described in 
iDuUemond fc DominikI (|2005l ). It considers growth un- 
der the extreme condition that only the particle of interest 
is settling toward the disc midplane and that all other 
particles remain small and suspended within the disc gas. 
In this limit 



dt 



^settle 



Zpg {Z) it^scttlc] (Tc 



(Bl) 



where mi is the mass of the particle, Z is the assumed dust- 
to-gas mass ratio within the disc, set to Z — 0.01 for all runs 
of our grain-growth models, pg (z) is the gas density local to 
height z above the midplane, Ca is the sound speed of the 
gas, and CTc is the coUisional cross-section of the growing 
particle, in this model set to CTc = t^'s'^i, where sai is the 
particle radius. The settling velocity of the grain, Usottio, is 
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Figure Bl. Depictions of grain growth in two cartoon models. 
The modeled, growing grain is highlighted in red, while the back- 
ground grains from which it grows are shown in blue with black 
outlines. Arrows depict the relative magnitude and direction (ran- 
dom in the random- motions model) of grain motions. Growth in 
the raindrop model is dependent on grain height, ^, above the disc 
midplane, while the grain growing in the random-motions model 
is assumed to exist always at 2: = 0. 
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radii. Because the inner regions of a disc are hotter and 
denser, particles there grow fastest and to the largest pos- 
sible sizes. Particles of a given size are initiated within 
the particle-transport simulations as soon as either of these 
grain-growth models grows a particle to that size. 



Figure B2. Grain growth using the two cartoon-grain-growth 
models. Blue: the raindrop growth model. Red: the random- 
motions growth model, p^ = 1 g cm^^. Growth here is modeled 
within the evolving fiducial disc model. 



defined for the case of Epstein drag. While the final, z = Q 
size of a grain in this model depends upon the height it is 
initiated, growth saturates for zq ~ 4_ffg. In our simulations, 
we therefore initiate grains from this height for the raindrop- 
growth calculation. 

For the random-motions growth model, we retain 
the assumption that all background grains remain small 
(sd2 = 0.1/J,m), allowing only the particle-of-interest to grow. 
Growth in this regime is controlled by the thermal velocity 
dispersion, Atjb, between grains, and is prescribed by 

— -— = Zpg,midAi;Bcrc , (B2) 

at 

Ai;b = 



8fcBT(mi +7712) 



7rmim2 

where 7712 is the mass of a background particle, and the 
growth model is run at the disc midplane where densities 
are highest and growth rates will be fastest. For this model 
we use the formal definition for the coUisional cross-section: 

ffc = TT (Sdl -I- Sd2) ■ 

Note that formally, growth via thermal velocity disper- 
sion is referred to as Brownian grain growth. However, in 
growth via Brownian motion, the entire small-grain popu- 
lation grows in size concurrently, so that growth is initially 
fast, but slows quickly as the thermal velocity drops respec- 
tive to the larger grain sizes. However, in turbulent models 
of grain-growth, growth via random motions can remain fast 
as the turbulent velocity dispersion increases for the larger 
size grains. See, e.g., the static-disc test cases p r esent ed in 
the appendix of iBirnstiel. Dullemond fc Braueij (|201Cl V By 
fixing the background grain population to small sizes in our 
random-motions grain-growth model, we retain some of the 
faster grain-growth experienced by larger grains in a tur- 
bulent environment, albeit, not to the extent that it occurs 
in formal turbulent-growth simulations and independent of 
the specific level of turbulence used within our disc-evolution 
simulations. 

In Figure IB2I we plot raindrop growth and random- 
motions growth within our fiducial disc model at three disc 
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